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Functions expressing dispersal probability decays with increasing distance are 
widely used in studies of animal movements. The inverse power function (IPF) 
exhibits the property of self-similarity, and hence should perform robustly 
against variation in marking efforts, allowing comparisons across studies. We in- 
vestigated this function property using dispersal data of Nickerl’s fritillary 
(Melitaea aurelia), a little studied checkerspot butterfly which is currently ex- 
panding in Central Europe. During mark-recapture in South Moravia, Czech Re- 
public, a single researcher worked for the entire flight period in 2005, while in 
2006 five researchers worked for just 5 days. Slopes of the fitted functions did not 
differ between the two seasons, illustrating the robustness of the function and 
suggesting the possibility to obtain reliable dispersal estimates even from brief 
marking campaigns. For both years, it was predicted that approximately one indi- 
vidual per one thousand would cross 10 km distance, the maximum distance sep- 
arating the most isolated colonies in the region. 
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1. Introduction 


The awareness of insect losses from habitat frag- 
ments surrounded by human dominated land- 
scapes, such as the intensive farmlands of West- 
ern and Central Europe, has raised interest in in- 
sect dispersal, or mobility related to gene flow 
(Habel et al. 2009, Hanski & Mononen 2011). 
Losses from habitat fragments are explicable by 
the metapopulation theory, which predicts that 
small local populations are prone to eventual ex- 
tinctions due to deterministic or stochastic pro- 


cesses, unless saved by occasional colonisation 
(Hanski 1999). Studies on model groups such as 
butterflies have provided ample evidence that in 
such isolated populations losses indeed occur 
(Wenzel et al. 2006, Leidner & Haddad 2011). 
The data amassed in empirical butterfly dispersal 
studies, using a variety ofapproaches from classi- 
cal mark-recapture to molecular and electronic 
tools, are now allowing broad generalisations 
which link dispersal to species’ life histories and 
aim to assist practical conservation (Hovestadt & 
Nieminen 2009, Stevens etal. 2010, Sekar 2012). 
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Checkerspot butterflies (Nymphalidae: Meli- 
teaeinae) represent a popular model in population 
ecology (e.g., McLaughlin et al. 2002, Hellmann 
etal. 2003, Ehrlich & Hanski 2004, Bulman et al. 
2007, Zimmermann et al. 2011). They often form 
spatially distinct colonies, are easily detectable 
both as larvae and adults, are rather poor fliers 
and occur in high densities. At the same time, they 
display aremarkable diversity of life history traits 
(Ehrlich & Hanski 2004), which transfers into 
varying levels of threat from intensive land use. 
Some European species are among the most en- 
dangered butterflies on the continent (e.g., Eu- 
phydryas maturna, cf. Freese et al. 2007), while 
others vary in status regionally (e.g. Melitaea 
athalia, cf. Hodgson et al. 2009). An interesting 
case in Central Europe is Melitaea aurelia 
Nickerl, 1850. This xeric grasslands specialist 
has declined severely during the 20th century 
(e.g., Van Swaay & Warren 1999, Beneš et al. 
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Fig. 1. Schematic rep- 
resentation of Zdanicky 
les hills, Czech Repub- 
lic, showing the study 
site, the distribution of 
xeric grasslands, and 
their occupancy by 
Melitatea aurelia within 
the mostly intensively 
farmed, or forested, 
area. 


2002), but is currently re-expanding to histori- 
cally occupied regions (Germany: Eichel & 
Fartmann 2008, the Czech Republic: Spitzer & 
Beneš 2010). It hence offers an excellent compar- 
ison with much more studied, but recently declin- 
ing, checkerspots. 

Mobility studies based on capture-mark-re- 
capture often miss rare long-range movements 
(Baguette & Schtickzelle 2003, Zimmermann et 
al. 2011), instrumental in the recolonisation/ex- 
pansion processes. Franzen and Nilsson (2007) 
recommended that a sampling area as large as 100 
km’ is necessary to estimate lepidopteran 
(re)colonisation ability. A way to overcome this 
logistic limitation is to use dispersal kernel func- 
tions, predicting long-range movements from 
short-range ones. This does not fully eliminate 
the problem with undetected long-range dis- 
persal, but at least allows comparisons among 
species, sexes, regions and seasons (Fric et al. 
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2010). Contributing to the debate (e.g., Baguette 
et al. 2000, Kuras et al. 2003) on the appropriate 
shape of dispersal kernel functions, Fric and 
Konvička (2007) showed that a particular func- 
tion form, the inverse power function (IPF), is 
particularly robust to variation in the marking ef- 
fort, owing to the scale invariance of the function 
slope; i.e., the parameter expressing decay of dis- 
persal with distance. 

The scale invariance of IPF offers an attrac- 
tive possibility to collect reliable dispersal data 
with reduced marking effort and to compare dis- 
persal data collected under varying study circum- 
stances. We examined this by performing a two- 
seasons mark-recapture study targeting move- 
ment patterns in a strong M. aurelia population 
inhabiting xeric grasslands in the Czech Repub- 
lic. In the first season, we followed a traditional 
approach, marking the population during the en- 
tire flight period. In the second season, we sub- 
stantially increased the study area, restricted 
marking to peak adult flight, and increased per- 
day marking intensity. Here, we (1) describe the 
basic demography and movement patterns in the 
population of M. aurelia, and (2) assess the ro- 
bustness of the movement estimates obtained 
while following different marking protocols. 


2. Material and methods 
2.1. Study system 


The study was carried out in a system of valleys at 
the southern slopes of the Zdanicky les hills 
(49°05’N, 17°02’E, altitude 245-330 m), South- 
em Moravia, SE Czech Republic. While most of 
Southern Moravia consists of intensively farmed 
lowlands, Zdanicky les is a ridge (maximum alti- 
tude 437 m), formed by base-rich Carpathian 
flysch and covered by oak-hornbeam and beech 
forests. Its south-facing valleys, historically cov- 
ered by a patchwork of gardens, hay meadows, 
orchards and pastures, were largely abandoned 
several decades ago, and now are in varying 
stages of succession. The prevailing vegetation is 
classified as Bromion erecti dry grasslands, Ge- 
ranion sanguinei forb-rich mantles, and Berbe- 
ridion xeric scrub (Chytrý et al. 2001). Despite 
high biotic diversity (e.g., 74 butterfly species re- 
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corded to date, including 17 nationally threatened 
ones, according to unpublished Czech Butterfly 
and Moths Recording data), complex land-hold- 
ing patterns have so far prevented legal protec- 
tion. A part of the dry grassland area is actively 
managed by conservation volunteers. 


2.2. Study species 


Melitaea aurelia is a univoltine butterfly, dis- 
playing colonial distribution at short-bladed xeric 
grasslands with a high cover of its host plant, 
Plantago media L. (Eichel & Fartmann 2008). 
Females lay small egg batches on the underside of 
host plant leaves, larvae are gregarious until au- 
tumn. Its total distribution stretches from France 
through Central and Eastern Europe, Transcau- 
casia, the Southern Urals and Siberia up to the 
Tian Shan region (Tolman & Lewington 2008). 
In the Czech Republic, a decline during the 20th 
century has restricted its distribution to the warm- 
est regions of the country, but recent records doc- 
ument a return of the species to parts of its former 
range. 


2.3. Mark-recapture data 


We worked in four adjoining valleys with a total 
dry grasslands area of 51 ha, divided into eleven 
subsites for the purpose of dispersal analysis (Fig. 
1). Two adjoining valleys (32 ha) were covered 
by marking in 2005, and two others (19 ha) were 
added in 2006. 

In 2005, the marking covered the entire adult 
flight period (25. VI-28. VII), every day with fa- 
vourable weather, from 10:00 to 17:00 (CEST). 
Each netted butterfly was assigned a numeric 
code and released at its point of capture. A single 
person carried out the marking, working for 32 
person-days in total. 

In 2006, the marking was limited to the peak 
of adult flight season (5.9. VII), but covered an 
expanded area. Five persons did the marking each 
day, resulting in the total marking effort of 25 per- 
son-days. Daily marking intensity, i.e. the num- 
ber of marking persons per day and unit area, was 
3.6 times higher in 2006 than in 2005 (5-fold in- 
crease in persons marking, but covering an area 
1.37 times larger). 
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2.4. Demography estimates 
and dispersal analysis 


Population size in 2005 was estimated using the 
Jolly-Seber method, POPAN parametrisation 
suitable for open populations with births, deaths, 
and migration, as implemented in the program 
MARK v. 4.3. The residence rate (Phi) combin- 
ing deaths and emigration, capture probability (p) 
and proportional recruitment (pent) are POPAN 
primary parameters, which may be constant for 
sexes and time (-) sex-dependent (g); or factori- 
ally (), linearly (T) or polynomically (Tx) time- 
dependent; additive (g + f) or interactive (g x f). 
POPAN derived parameters are daily recruitment 
(Bi), daily population size (Ni), and total popula- 
tion size (Ng) (Schtickzelle et al. 2002, Vlasanek 
et al. 2009). 

To compare mobility between years and 
sexes, we used data on patch-to-patch move- 
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Fig. 2. Daily numbers of 
captures of Melitaea 
aurelia males and fe- 
males. — a. In 2005. 

— b. In 2006. The ar- 
rows in panel (a) show 
the 2005 days pheno- 
logically corresponding 
to the five marking days 
in 2005. 


ments, lumped them into 50-m distance classes 
(e.g., Drag et al. 2011), and fitted the inverse 
power function, expressing the probability den- 
sity (J) of movements as 


I= C.D” (la) 
or, in logarithmic form, 
lIn Z= In C- n(ln D) (1b) 


To compare the 2006 (shorter term, expanded 
area) data with the 2005 data, we used two 2005 
datasets. One consisted of all available data. For 
the second, we selected movements on five con- 
secutive days with a sex ratio corresponding to 
2006 data (see Fig. 2a). 

We fitted IPFs on these three datasets and 
compared obtained slopes and intercepts using 
the Tukey HSD test (Zar 1996). 
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Table 1. Overview of mark-recapture data collected during the study for Melitaea aurelia males (M) and females 


(F). 











Marking Individuals Individuals Capture Movements 
days M/F captured M/F recaptured M/F events M/F M/F 
2005 32/29 304 / 108 318 / 98 649 / 230 146/20 
2005 reduced 5/5 113/55 62/11 191/67 29/8 
2006 5/5 827 / 468 455 / 110 1375 / 598 274/151 
3. Results 1). The total population size for 2006 was not re- 


For the 2005 marking (Table 1), the best demog- 
raphy model was Phi(t)p(g+t)pent(g*72)Ng), 
containing 50 parameters (AICc = 2162.58). The 
second-ranking candidate model differed by 
AAICc > 3.0. Melitaea aurelia residence Phi thus 
responded factorially to the day of marking, inde- 
pendently of sexes; capture probability p varied 
among marking days and was always higher in 
males; and recruitment pent was approximated by 
domed quadratic polynomials, differing in shape 
between sexes and revealing protandrous recruit- 
ment (Fig. 3). The total estimated population size 
was 490 (+35 SE) males and 270 (+36 SE) fe- 
males. 

In the 2006 marking, we marked 3.1 times 
more individuals than in 2005. Based on sex ratio 
on capture, the marking period in 2006 (5-9. 








alistically identifiable from the data, due to too 
low number of marking days, but probably 
reached a few thousand individuals. 

The mean / median flight distances in 2005 
were 307436 SE m/ 200 m for males and 203+6 
m/ 100 m for females; in 2006, 377423 m/200 m 
for males and 398+8 m/250 m for females. Based 
on the reduced 5-day data, the mean / median dis- 
tances were 256+10 m / 200 m for males and 
13244 m/ 100 m for females. In both years, the 
longest detected movements were those of males 
(Table 2). 

All fitted IPFs approximated the dispersal 
data significantly (Table 2), with all R’ > 0.7 (Fig. 
4). Comparing the functions across years and 
sexes did not reveal any overall difference in in- 
tercepts (F = 0.25, df= 1, 5, p = 0.618) or slopes 


























VII) corresponded to 10.-14.VII in 2005. The 
males had already passed their annual peak, 
whereas females were just peaking (Fig. 2). The 
number of butterflies marked in those phenologi- 
cally corresponding 2005 days was 168 (Table 


Estimated population size 


(F=0.08, df= 1,5, p=0.775). In other words, the 
short-term marking in 2006 detected identical 
dispersal parameters as the entire flight period 
marking in 2005. Post-hoc comparison of all pos- 
sible combinations, however, suggested a differ- 
ence in slopes between the total data of females in 
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Table 2. Test statistics for the inverse power functions linear regressions and derived estimates of Melitaea 
aurelia movement probabilities to distances longer than those covered by the marking. 











d.f.; F; p Predicted movements (%) 
Max. 
1km 2km 5 km 10 km distance 
2005 
males 1,15; 66.4; <0.001 2.94 1.03 0.25 0.09 1,250 m 
females 1,7; 38.7; 0.004 3.12 1.36 0.46 0.20 600 m 
2005 reduced 
males 1,11; 20.8; <0.001 2.49 1.14 0.41 0.18 750 m 
females 1,3; 47.8; 0.002 1.06 0.52 0.2 0.01 450 m 
2006 
males 1,13; 51.4; <0.001 1.79 0.73 0.22 0.09 2,000 m 
females 1,14; 38.2; <0.001 4.18 2.26 0.99 0.54 1,050 m 





2005 and males in 2006 (at p = 0.05), but this sin- 
gle difference was likely influenced by the differ- 
ent sample sizes in the two regressions. 

The IPF-derived predictions of movements to 
5 (10) km distances suggest that a population con- 
taining 1,000 (10,000) individuals will likely 
contain a few individuals moving to such dis- 
tances. The occupancy pattern in a wider region 
(Fig. 1) shows that patch-to-patch distances 
among occupied xeric grassland patches are with- 
in this distance range. 


4. Discussion 


Estimation of the parameters of the inverse power 
function dispersal kernel using the M. aurelia 
mobility data obtained in two consecutive years 
resulted in very close function parameters, and 
hence predicted identical migration rates. The 
functions were identical despite major differ- 
ences in data collection in the two years: while in 
2005 we covered the entire flight period, in 2006 
we marked the butterflies only during the peak 
flight, albeit with higher marking intensity. 

The unchanged function performance sup- 
ports the conjecture by Fric and Konvicka (2007) 
that given the scale invariance of IPF, it should 
predict identical dispersal rates under a great vari- 
ety of marking protocols. The authors showed 
this by subsampling data for two species, Euphy- 
dryas aurinia (Rottenburg, 1775) and Parnassius 
mnemosyne (Linnaeus, 1758). Here, we in fact 
subsampled the 32-day 2005 data to 5-day 2005 


data. On the other hand, the increased marking in- 
tensity in 2006 expanded the data set, but none of 
the two procedures changed the original (i.e., to- 
tal data 2005) IPF kernel. 

Melitaea aurelia adult demography patterns 
do not deviate from those found in other checker- 
spots (Ehrlich & Hanski 2004, Konvička et al. 
2005, Schtickzelle et al. 2005, Hodgson et al. 
2009), or other single-brood nymphalids (e.g., 
Proclossiana eunomia: Schtickzelle et al. 2002; 
Boloria aquilonaris: Baguette & Schtickzelle 
2003). Neither does the increase in numbers be- 
tween 2005 and 2006 appear unusual, as abun- 
dance fluctuations by an order of magnitude oc- 
cur frequently in butterfly populations (Hellmann 
et al. 2003, Schtickzelle et al. 2005, Cizek & 
Konvička 2009). In species with communal lar- 
val nests, such changes may be particularly prom- 
inent as entire broods are easily destroyed by ca- 
tastrophes, parasitoids or pathogens (Ehrlich & 
Hanski 2004, Schtickzelle et al. 2005), which 
might be responsible for the rapid patch occu- 
pancy turnover in checkerspot metapopulations 
(Hanski 1999). 

Mean and median movements also fell within 
the ranges reported in previous checkerspot stud- 
ies i.e. a few hundred meters for mean distances 
and 1-2 kilometres for maxima (e.g., Wahlberg et 
al. 2002, Baguette 2003, Hovestadt & Nieminen 
2009, Zimmermann et al. 2011). The 2006 mark- 
ing returned longer mean distances, which could 
have been caused by either the increase of mark- 
ing area (e.g., Schneider 2003, Zimmermann et 
al. 2011), or by density-dependent emigration 
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Fig. 4. Fitting of inverse power functions (IPF) expressing the probability density of Melitaea aurelia movements 
to certain or longer distances against the distances (/n-In scale). Short-term intensive marking in 2006 returned 
identical dispersal kernels (i.e., slopes and intercepts of the function) as did the marking in the whole flight period 


in 2005. 


(Nowicki & Vrabec 2011); the data do not allow 
distinguishing between the two alternatives. The 
mean flight distances obtained with the reduced 
2005 data were by one-third (males) or even by 
one-half (females) shorter compared to the total 
2005 data, which is attributable to a lower chance 
of detecting long individual displacement during 
the shorter marking time. It is hence more encour- 
aging that IPF did not differ among the three data 
sets. 

The ranges of fitted IPF slopes (—1.1 to —1.5 
for males, —0.9 to—1.2 for females) were also sim- 
ilar to other checkerspot studies (e.g., Fric et al. 
2010, Zimmermann et al. 2011). The overlap of 


the values with other checkerspots suggests a tax- 
onomic conservatism in mobility. This justifies 
using data from closely related species, such as 
various checkerspots, for such purposes as mod- 
elling or reserve design (cf. Wahlberg et al. 1996, 
Fric et al. 2010), however, the further verification 
on other taxonomic groups is needed. Differences 
between sexes were not statistically significant, 
although male slopes tended to be steeper, indi- 
cating fewer long-distance movements compared 
to females. This pattern, explicable by different 
fitness consequences of dispersal for sexes 
(Zonneveld & Metz 1991), is frequent in butter- 
flies (Scott 1975, Baguette & Néve 1994, Hanski 
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& Thomas 1994, Gutiérrez et al. 1999). Female 
dispersal results in establishment of new colo- 
nies, and ultimately, in range expansions. Identi- 
cal movements in butterfly sexes, however, were 
also reported for some checkerspots (Wahlberg et 
al. 2002, Wang et al. 2004). 

In Germany, M. aurelia often forms small 
colonies at calcareous grassland islets (Eichel & 
Fartmann 2008). Habel et al. (2009) showed that 
despite living in small colonies, it maintains high 
levels of genetic polymorphism by occasional 
long-distance dispersal. Our IPF-derived move- 
ment probabilities (Table 2) predict that one indi- 
vidual per thousand reaches a 10 km distance. 
Given the xeric grasslands distribution in the re- 
gion (Fig. 1), this should maintain a mutual inter- 
connection of the system. Recent capture-mark- 
recapture studies at scales of hundreds of square 
kilometers (Franzen & Nilsson 2007, Zimmer- 
mann et al. 2011) revealed that occasional long- 
distance dispersal movements occur even in ap- 
parently sedentary species. Such movements are 
obviously rare, but critical for species survival in 
increasingly fragmented landscapes. 

In summary, this study supports the earlier 
observation by Fric and Konvička (2007) that 
dispersal probabilities obtained with the inverse 
power function should be robust with respect to 
variation in marking protocols. While the earlier 
study re-sampled single-year data, this study ana- 
lysed data sets generated under different marking 
protocols. Despite unavoidable year-to-year vari- 
ation in such factors as weather (effects on butter- 
fly mobility: e.g., Cormont et al. 2011) or popula- 
tion density (e.g., Nowicki & Vrabec 2011, 
Konvička et al. 2012), the dispersal probability 
decays did not change. This opens the exciting 
possibility of obtaining reliable dispersal esti- 
mates even from brief marking campaigns. Al- 
though we understand that a comparison based on 
just two seasons should not allow too much opti- 
mism, the possibility still deserves investigation. 
We also understand that too brief marking cam- 
paigns cannot provide information on such popu- 
lation traits as longevity, or temporal develop- 
ment of sex ratio. Regarding dispersal, however, 
it is increasingly apparent that to predict species’ 
fates in fragmented landscapes, the knowledge of 
habitat patch distribution and occupancy patterns 
may be more critical than movement data (cf. Fric 
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et al. 2010). Dispersal is a probabilistic process 
and the chance to reach more distant sites in- 
creases with abundance at the source colonies, 
highlighting the need for appropriate conserva- 
tion management of as large as possible stretches 
of suitable habitats (Eichel & Fartmann 2008). 
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